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Pure gauge compact QED on hypercubic lattices is considered with periodically closed monopole currents 
suppressed. We compute observables on sublattices which are nested around the centre of the lattice in order 
to locate regions where translation symmetry is approximately recovered. Our Monte Carlo simulations on 24*- 
lattices give indications for a first-order nature of the U(l) phase transition. 
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1. INTRODUCTION 

The presence of metastabilities at the phase 
transition of compact QED became obvious in 
simulations on 6*-lattices with periodic bound- 
ary conditions (p.b.c.)||l|. With increasing lat- 
tice size, there was growing confidence that the 
U(l) phase transition is exhibiting a weakly first- 
order nature, culminating in the observation of 
extremely strong hysteresis effects appearing in 
Monte Carlo series on 16''-lattices|Q]. 

This widely accepted opinion has been ques- 
tioned recently. The authors of Refs. have 
implemented compact U(l) on a lattice with triv- 
ial homotopy group, i.e. SH^, the surface of H^. 
No indications for metastabilities have been found 
at f3c in simulations on lattices of different sizes — 
up to a number of degrees of freedom equiva- 
lent to a 16^-lattice with hyper-toroidal topology. 
Similarly, fixed boundary conditions (f.b.c.) hint 
at a continuous transition ||^,^. 

The implementation of U(l) on SH4 as well 
as with f.b.c. explicitly breaks the translation 
symmetry and produces inhomogenities of the 
free energy. One might argue that the absence 
of first-order effects can be explained in these 
terms. Alternatively it has been suggested in 
Ref. 0] that the existence of topologically non- 
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trivial monopole loops wrapping around the lat- 
tice due to p.b.c. might lead to a first-order sig- 
nature at /3c- 

We can check this assumption by studying 
the model with p.b.c. but periodically closed 
monopole loops forbidden and looking for first- 
order signals. For this purpose, we suppress 
monopole currents emerging in the Monte Carlo 
update that would point outside the boundary 
shell and thus would be closed due to p.b.c. 

As monitoring quantity we measure the average 
plaquette on the hyper-surface of subcubes nested 
around the centre of the lattice and we try to 
establish a region where translation symmetry is 
recovered approximately. 

2. SUPPRESSION OF MONOPOLES ON 
BOUNDARY 

We suppress monopoles on the boundaries in 
such a way that spinwaves can penetrate them. 
The monopole current is defined byj^ 



npa(x + i/) - np„{x) 
2,-1,0,1,2, 



(1) 



where the number of Dirac sheets n^i, [x) is given 
by the difference 
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Figure 1. Plaquettes 
contributing to m^(x). 



Figure 2. Measure- 
ments on nested shells. 



between physical flux $,^1/(2;) and plaquette flux 
For a link Up, to be updated in a 
boundary-cube with normalvector e^, we search 
for the plaquettes associated with the Dirac 
sheets npa{x) that contribute to the monopole 
current m^(x) being parallel to e^. These plaque- 
ttes are depicted in Fig. |l| with dotted lines. If, 
for any of the corresponding cubes, the change in 
flux is larger than 27r, the proposed link change 
is rejected in the Monte Carlo. This particular 
update rule fulfills detailed balance and is equiv- 
alent to an infinitely large chemical potential for 
the monopoles on the boundary shell||^ with the 
action S — Sw + S\ . Sw is the standard Wil- 
son action and S\ is the additional term which 
suppresses monopole currents at the lattice 
slice Xp = 0: 

Sx = X-^ ^ I TO^(x) I, A ^00. (3) 

S\ explicitly breaks translation symmetry. 

3. RESULTS 

On the 16^- as well as on the 24''-lattice, we first 
scanned the /3-range in order to get a rough idea 
about the location of the (pseudo)-critical region. 
It is shifted towards smaller /3-values as expected 
from monopole suppression. Table 1 shows the 
number of Metropolis sweeps we performed sub- 
sequently for the such chosen /3-values and the 
plaquette averaged over the whole lattice. On 
the 24'*-lattice in the average about 1% of the 
proposed updates are rejected at the boundary. 

The plaquettes have been sampled on the sur- 
faces of sub-hypercubes nested around the centre 
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Figure 3. Average plaquette (dotted line) and 
plaquettes on shells 1,3,5 and 7 (from top to bot- 
tom) on the 16^-lattice. 

of the lattice. The 3D projection of these shells is 
shown in Fig. ^. The boundary shell is denoted 
by the index 1, the inner shell gets the index L/2. 

Table 1: Lattice sizes, /3-values and # of sweeps. 



size 


/3 


# of sweeps 


< P > 


16^ 


0.990 


30000 


0.6290(2) 




0.995 


10000 


0.6408(4) 




1.000 


60000 


0.6530(2) 




1.002 


60000 


0.6584(1) 




1.003 


60000 


0.6599(1) 




1.005 


60000 


0.6630(1) 




1.010 


10000 


0.6688(2) 




1.0060 


60000 


0.6546(1) 




1.0065 


60000 


0.6572(1) 




1.0070 


250000 


0.65856(5) 




1.007125 


200000 


0.65935(9) 




1.0075 


80000 


0.66010(9) 




1.0078 


50000 


0.66099(6) 




1.0080 


55000 


0.66139(5) 



The boundary shell is always in a deconfined 
state. This is a well known phenomenon in the 
absence of monopoles Therefore, the varia- 
tion of the plaquette value within the mctastable 
region is very small for this shell. As the bound- 
ary shell contains the outweighting part of the 
sites, the plaquette average over the whole lattice 
is dominated by the contribution of the boundary. 
We found the influence of the boundary decreas- 
ing towards the centre, i.e. the plaquette values 
show increasing 'steps' (as functions of /3) going 
to the inner shells, see Fig. || for the 16^-lattice. 
For inner shells the plaquette values are close to 



3 



shell 1 




Figure 4. Time series on 24'^-lattice for different 
shells. 

each other and their behaviour is more and more 
like that of the plaquette of the H4 system with 
p.b.c. We interpret this results as the onset of 
the recovery of translation invariance around the 
centre. 

We carefully analyzed the time series for the in- 
ner shells on the 16^-lattice but the series merely 
exhibited a chaotic behaviour and no metastabil- 
ities could be identified. Therefore we decided 
to work on a very large 24^-system and we per- 
formed altogether nearly one million sweeps. 

On this lattice, we found a much broader region 
in 4-d space at the centre where the behaviour of 
the system resembled that of U(l) on H4 with 
p.b.c. and we found evidence for metastabilities 
in the time series. In Fig. ^, the time series of 
plaquettes for the shells 1, 5 and 9 are depicted 
for /? — 1.007125. At this value of /3 we encoun- 
tered long metastabilities extending over more 
than 10000 sweeps as can be most clearly seen 
in the series of shell 9. As the fluctuations of the 
plaquette in inner shells are larger than the gap 
between the two metastable states we performed 
a windowing procedure in order to average out 
the fluctuations on small scales and applied the 
method to the time series of shell 9. The result is 
presented in Fig. ^, where we see a non-ambiguous 
double-peak structure of the histogram. More de- 
tailed results will be presented in a forthcoming 
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Figure 5. Histogram of shell 9 after removing of 
fluctuations. 

paper. 

4. CONCLUSIONS 

In our simulations 

- we did not observe flrst order signals on small 
lattices because of boundary effects. We believe 
this is true for f.b.c. as well. 

- we observed metastabilities coming up on the 
24^-lattice. This implies that periodically closed 
monopole loops can not be the reason for the flrst 
order nature of U(l) with p.b.c. 
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